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Abstract 

We develop a framework for the reconstruction of the non-forward kernels which govern the 
evolution of twist-two distribution amplitudes and off-forward parton distributions beyond leading 
order. It is based on the knowledge of the special conformal symmetry breaking part induced by 
the one-loop anomaly and conformal terms generated by forward next-to-leading order splitting 
functions, and thus avoids an explicit two-loop calculation. We demonstrate the formalism by 
applying it to the chiral odd and flavour singlet parity odd sectors. 
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1 Introduction 



The possibility to access new characteristics of hadrons by means of the deeply virtual Compton 
scattering |2], |3| and the hard diffractive hadron electroproduction [[| |J processes has recently 
initiated a growing phenomenological interest in the underlying non-perturbative elements — the 
so-called off-forward parton distributions (OFPD) — which parametrize hadronic structure in 
these reaction making use of the QCD factorization theorems. The main feature of the processes 
is a non-zero skewedness, i.e. plus component, A + = rj, of the t-channel momentum transfer A. 

One of the central issues which has been addressed in this context is the description of the 
scaling violation phenomena in the cross section via the evolution of the off-forward parton dis- 
tributions. Since the OFPD is defined as an expectation value of a non-local string operator, its 
Q 2 -dependence is governed by the renormalization of this operator Fourier transformed to the 
momentum fraction space. Inasmuch as the generalized skewed kinematics can be unambiguously 
restored || from the conventional exclusive one, known as Efremov-Radyushkin-Brodsky-Lepage 
(ER-BL) region rj = 1, in what follows we deal formally with renormalization of the ordinary 
distribution amplitudes which obey the ER-BL equation ||. [?J 

—^<t>(x, Q) = V (x, y\a s (Q)) I <f>(y, Q). (1) 



dlnQ 

Here we have introduced the exclusive convolution 

f 1 

® = / dy, 
Jo 

to distinguish it from the inclusive one used later. Here <fi = C^jj) is the two-dimensional vector 
and V(x,y\a s ) is a 2 x 2-matrix of evolution kernels given by a series in the coupling. 

Several methods have been offered so far to solve the off-forward evolution equation: numerical 



integration ||, expansion of OFPD w.r.t. an appropriate basis of polynomials || [K|, mapping to 
the forward case0 |12|, [13] and solution in the configuration space [14], |I5| . The last three methods 
are based on the well-known fact that operators with definite conformal spin do not mix in the 
one-loop approximation. Beyond leading order the latter two methods can only be applied in the 
formal conformal limit of QCD where the /3-f unction is set equal to zero and making use of the 
conformal subtraction scheme which removes the special conformal symmetry breaking anomaly 
appearing in the minimal subtraction scheme. Thus, we are only left with the former two methods 
which allows for a successive improvement of the perturbative approximations involved. 

Up to now only orthogonal polynomial reconstruction method has allowed the analysis of 
the scaling violation in the singlet sector in two-loop approximation since only the anomalous 



2 This idea has earlier been applied directly to the kernels in 



1 



dimensions required in the formalism were available so far [|TT| , [T6| , p~T|] . This was sufficient to get a 
first insight into the NLO evolution corrections. However, in order to have an access to the whole 
kinematical region, especially for small x, rj and high precision handling of the x ~ rj domain, one 
should look for a more efficient numerical treatment. This can be achieved with the first method 
alluded to above. To do the direct numerical integration of the evolution equation one needs 
the corresponding evolution kernels whose Gegenbauer moments define the anomalous dimensions 
mentioned earlier. For the time being the former were available at LO order only. The flavour 
non-singlet ER-BL kernel (rj = 1) was obtained in NLO by a cumbersome analytical calculation 
15] , |19t P0| . As we have mentioned above the continuation to r] G [0, 1] is a unique procedure 0, 



so that one can obtain in a simple way the evolution kernels for OFPD. The goal of this paper is 
to outline a method that allows one to construct the singlet ER-BL kernels by applying conformal 
and supersymmetric constraints where the latter ones arise from the J\f = 1 super Yang-Mills 



theory pT| , |22|| . In this way we can avoid the direct diagrammatical calculation which would be 
very difficult to handle otherwise since no appropriate technology has been developed yet. 

The derivation is based on the fairly well established structure of the ER-BL kernel in NLO. 
Up to two-loop order we have 

2 



V(x,y\a s ) = %LV®{x,y) + V^(x,y)+0(a 



2tt 



a, 



2n 



(2) 



with the purely diagonal LO kernel in the basis of Gegenbauer polynomials and NLO one 
separated in two parts: V^\x, y) = V D( - 1 \x, y) + V ND ^(a;, y), with the diagonal part which is 



entirely determined by the well-known forward DGLAP splitting functions P(z) [|TT 



! V (x,y) 



o Ni(y) Cj 



(2x-l)z jAB P(z)C" {B) {2y-l) 



(3) 



where NAv) 



-4^+i r 2 (^)r(2^+i) 



and w(y\v) = (yy) u are the normalization and weight factors, 



respectively. The non-diagonal piece is fixed completely by the conformal constraints |T7j 

A) 



V^\x,y) = -(l-V) V. 



y(0) + ^ j 
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v(°) 



fay), 



(4) 



in terms of 



C F WVW 2T F N f Q G V® 
C F GQ V W C A GG V^ 
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C F QQ g 
r< gq„ r< gg„ 



(5) 



the ER-BL kernels at LO and the special conformal symmetry breaking matrix g. Here (3q = 
^TpNf — yC^ is the first expansion coefficient of the QCD ^-function. In the parity odd sector 
the dotted kernel, V, is simply given by a logarithmic modification of the . Due to subtleties, 
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appearing in the parity even case |TT]] , we deal here, for the sake of simplicity, only with the parity 
odd and transversity sectors. 

The main problem is thus to restore the diagonal part of the NLO kernels. Since the use of 
Eq. (|3|) beyond LO is extremely complicated in practice, we are forced to look for other solutions. 
It turns out that the bulk of contributions in the ER-BL kernel can be deduced by going to the 
forward limit making use of the reduction 

P(z)=LIMV(x,y) = \im±\ I (-.-). (6) 




Then the difference^ 

P( Z ) - pcross-laddcr^ _ JJM V N °(x, V ) 

can be represented in terms of inclusive convolutions of simple splitting functions and the back 
transformation to the exclusive kinematics is trivial. The contributions of the purely diagonal 
cross-ladder diagrams y cross - ladder ( z ) can be found from the known QQ sector |1| [19|, |(J exploit- 
ing the Af = 1 supersymmetric constraints. 

The paper is organized as follows. In the next section we analyze the structure of the known 
flavour non-singlet ER-BL kernel and state the benchmarks of the formalism. The structure 
observed will give us a guideline for construction of all other kernels: quark chiral odd sector is 
considered in Section 3 and parity odd flavour singlet one is discussed in Section 4. Finally, we 
give our conclusions and an outlook . 



2 Structure of ER-BL kernel in non-singlet sector. 



It is very instructive to demonstrate the machinery in the simplest case of non-singlet sector. Since 
the explicit two-loop calculation is available |1| [19], [2(J] the direct comparison can be made. The 
NLO QQ-kenael can be decomposed in colour structures asQ 



V{x,y\a s ) 



a 



Z7T 



+ 



2ttJ 



C F V F (x, y) - y v p(x> V) - {Cf - ~y ) >! (•'(•''■ n) 



(7) 



3 Here V ND is understood without the (I — 2?)-projector. 
4 We omit the superscript QQ later in this section. 
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with the LO kernel V^°'(x,y) = [v(x,y)] + , where 

— x ( 1 

v(x,y) = 6(y - x)f(x,y) + 6(x - y)f(x,y), and f(x,y) = -[l + 



y \ y-xj 

The shorthand notations x = 1 — x and / = f(x, y) are used throughout the paper. The "+"- 
prescription is conventionally defined by 

[V{x, y)} + = V(x, y) - 5{x - y) j dz V(z, y). 

Let us now recall a few properties of the kernel that are useful for the following considerations. 
Due to absence of the conformal symmetry breaking counterterms at leading order for the renor- 
malization of the composite operators with total derivatives, one can use its consequences to fix 
the eigenfunctions which turn out to be the Gegenbauer polynomials C^ 2 (2x — 1) f| 0. Thus, 
the LO kernel is symmetric with respect to the weight function xx: yyV^°\x,y) = xxV^°'(y,x). 
Its eigenvalues are given by the anomalous dimensions appeared in the analysis of deep inelastic 
scattering. Thus, it is not surprising that a simple limit already mentioned in Eq. (^) gives us the 
DGLAP kernel §: 

P(z) = LIM V(x, y) = lim -^-V cxt (-, -) . (9) 



t^O x \x X 

To perform this limit, we have to extend at first the ER-BL kernel, originally defined in the domain 
< x, y < 1, to the whole region x, y G (— oo, oo) by a unique procedure which is given in practice 
by the replacement, e.g. at leading order, of the ^-function by 

9(y-x)^e(l-^y(^)sign(y). (10) 

If a kernel is diagonal in the ER-BL representation, we can restore it from the known DGLAP 
kernel by the integral transformation Because of branch cuts appearing in the convolutions 
of the NLO terms with the transformation kernel, it is highly nontrivial to handle the inverse 
reduction to the exclusive kinematics. 

At NLO the kernel (0) contains besides a pure diagonal part with respect to the Gegenbauer 
polynomials also a non-diagonal part located in Vp(x,y) and Vp(x, y). These parts are predicted 
by conformal constraints (see QQ-entry of Eq. @f) and are fixed by the one- loop special conformal 



X 



anomaly kernels |TT], [L6], [17| : 

v(x,y)=6(y-x)f(x,y)\n- + \ X ^ X ), g(x, y) = -9(y - x) ^— *l + \ X ^ *\ . (11) 

y [y^yj y-x [y^yj 

Let us now analyze in detail the contributions to the NLO kernel from different colour struc- 
tures. The expressions for C\ terms arise from Feynman diagrams containing quark self-energy 
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insertions and ladder graphs^. In order to subtract the ultraviolet (UV) divergences in subgraphs 
it requires the LO renormalization of the composite operator to which these lines are attached to. 



The explicit calculation gives [18, 19, 20 



V F (*,V) - %-x){(l-2C(2))/ + 3^-g/-|).„^-(/-7).n = l„(l-^ 

+ ( f + l^tf £ 1 _ ijnifl + In a; - 2 In x) + [ X ~* X 1 . 
\ 2yJ y\ 2y \y -> y I 



'121 



Making use of the known non-diagonal part (^), Vjr can be represented up to a pure diagonal 
term, denoted as D F (x,y), by the convolution 



V F {x,y) 



c c e \ 

v®t> + gCg>t> — t><g>g (x, y) + D F (x, y). 



(13) 



To find an appropriate representation of this missing diagonal element we first take the forward 
limit. Since the forward limit of the convolution isPl 



LIM [v}+ ®[v}+ = {LIM [«]+} ® {LIM [v]+} , 



(14) 



where we have introduced the inclusive convolution 

P 1 (z)®P 2 (z) = ( dx ( dy5(z-xy)P 1 (x)P 2 (y), 
Jo Jo 

and the commutator g ® — g drops out in the forward limit, we obtain 

UMV F (x,y) = -p®p + LIMD F (x,y), 



(15) 



where p = LIMv = p(z) \nz + 1 — z and p(z) = LlMv(x,y) = (1 + z 2 )/(l — z). The comparison 



of LIMV F (x,y) with the corresponding part of the DGLAP kernel [23 

4 



P F {z) = !^-- 2((2) - -\nz + \n 2 z -2\nz\n(l- z)~jp(z) 



1 — 3z . l + .z l2 
+ 1 — z H mz In z, 



(16) 



yields the result in which all double log terms are contained in the convolution p®p and, therefore, 
only single logs survive in D F (z) = LlMD F (x,y): 



DAz) 



i 

P F +p®p 

-^p a (-z) \nz- p a (z) | In ^ - 2 ln(l - z) - - 
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p(z) 



(17) 



5 For simplicity wc imply the diagrams in the light-cone gauge jl8], [H 
6 We remind as well that [A]+ ®[B] + = [C] + . 
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Here we have introduced for convenience the kernel p a (z) = 1 — z. The next important point is 
that the remaining log terms can be represented as convolutions of p a and p. Thus, we have finally 

D F (z) = ±p«®{2p + p a } (z) + ±p{z) + h -p a {z). (18) 

Since Dp(x,y) is by definition diagonal, the extension of Dp(z) towards the ER-BL kinematics is 
trivial: 

D F (z) ^D F {x,y) = l -v a ® (2v + v a ) (x,y) + ^v(x,y) + ^v a (x,y), (19) 

where a new diagonal element is v a (x,y) = 9{y — x)| + Qix — Evaluating the convolutions 
one can establish the equivalence of our prediction with Eq. (|12}). 

Next, the Feynman diagrams containing vertex and self-energy corrections provide Vp pro- 
portional to (3o- Its off-diagonal part is induced by the renormalization of the coupling and is 



contained in the dotted kernel (11 



V (x,y)=v(x,y) + Dp(x,y). (20) 
The remaining diagonal piece, Dp, is deduced from the known NLO DGLAP kernel 



Pf>{z) = ~p(z)+p a (z)+p(z) (21) 
by going to the forward kinematics and restoring then the missed contributions from it. Thus, 

Dp(x,y) = -v(x,y) + v a (x,y). (22) 



Indeed, the final result coincides with |T8|. [19]. g0| . 

Finally, we come to the contribution which mainly originates from the crossed ladder diagram 
proportional to (Cp — Ca/2)' 

V g (x, y) = 2v a (x, y) + ^v(x, y) + (g(x, y) + j ^ ^ . (23) 

Since this diagram has no UV divergent subgraph and thus requires no subtraction, its contribution 
has to be diagonal w.r.t. the Gegenbauer polynomials. This is obvious for the first two terms 
appearing in Eq. (p3|). The function^ G(x,y) contains in the unusual 9(y — x)-structure the 
mixing between quarks and antiquarks 

G(x, y) = 6(y - x)H(x, y) + % - x)H(x, y), (24) 



7 We have slightly changed the original definition given in |l8|, by G(x,y) + 26(y — x)f \ny\nx — > G(x,y). 



6 



with 



H(x,y) 
H(x,y) 



f(U 2 (x) + \ny\nx) - fLi 2 (y) 
(f - /) (Li 2 In 2 y)+f (Li 2 (y) - U 2 (x) - \ny\nx) 



(25) 
(26) 



where Li 2 is the dilogarithm. It can be easily checked that the G-contribution (not the terms H 
and H separately) is symmetrical w.r.t. the weight xx. Performing the limit (|9|) we obtain the 
following correspondence with the non-singlet DGLAP kernel PJ: 



where 



G(z) = LlMG(x,y) = 6(z)6(l - z)H(z) + B{-z)B{\ + z)H(z), 

H(z) = UMH(x,y) = p(z) (in 2 2 - 2((2)) +T(z), 
H(z) = UMH(x,y) = 2p(z)S 2 (-z) +T(-z). 



(27) 



(28) 
(29) 



Here S 2 (z) = f In ^ and T(z) = 2(1 + z) In z + 4(1 - z). We should emphasize that 

this structure of G is the most general, especially, it is the only contribution that contains Spence 
functions. In the forward limit we obtain therefore a typical combinations given in (p8[) and (p9|) , 
which can be found in all other channels as well. This observation provides us with a hint for the 
construction of all singlet G kernels in the ER-BL representation. For completeness, we give the 



corresponding part of the DGLAP kernel [£3 

P G (z) = UMV G (x,y) 



2p a (z) + -p(z)+G(z), 



(30) 



and G(z) defined above in Eqs. 

Recapitulating the results obtained in this section, we have observed a rather simple structure 
of the non-singlet NLO kernel in the QQ-channel. Up to the diagonal G-function, which is in 
fact the only new element in the two-loop approximation, we can represent all other terms by a 
simple convolution of LO kernels already known. It is not accidental but a mere consequence of 
the topology of contributing Feynman graphs at 0(a 2 s ). Thus, we anticipate the same feature to 
appear in all other channels as well. 



3 Quark kernel in chiral odd sector. 

After we have outlined and tested in the preceding section our formalism, we can apply it to the 
previously unknown transversity two-loop ER-BL kernel. We decompose the transversity kernel 
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analogous to the chiral even case ([?[). We also use the same decomposition for the DGLAP kernels 



24 1 - The leading order kernel is 



with 



V^ T (x,y) = [v b (x,y)] + - -5(x - y), 



v b {x,y) =9(y-x)f b (x,y) + 0(x-y)f (x,y), and f b (x,y) 



x 1 



(31) 



(32) 



yy-x 

The non-diagonal part has been analyzed in Ref. [T7| and is completely analogous to the chiral 
even case discussed above. Thus, 



® V^ T + [g] + & V^ T - V^ T & [g] + (x, y) + D T F (x, y) 



(33) 



where v b is obtained from Eq. flTT| ) by replacing / by f b and the g kernel is the same as in Eq. 
(|Tl"D. Taking the forward limit of V F (x, y) and comparing it with the known result for the DGLAP 



kernel in Ref. p4 |, we find the following trivial representation of the remaining diagonal part 



D T F {x,y) 



2 r 



v b (x,y) 



19 
24 



5{x - y). 



(34) 



There is essentially no extra work required to find the contribution proportional to the /Jo- 
function, since it can be easily traced from the DGLAP kernel to be 



5 r 



v /3 (x,y) = - v (x,y) + v (x,y) 



13 
12 



6(x - y). 



(35) 



The case of the G function is easy to handle as well. If we replace / by f b in the definition 
we obtain the diagonal G T (x,y) kernel. Taking the forward limit and comparing it with 
the DGLAP kernel, we immediately find the remaining ^-function contribution, so that the whole 
result reads 



V£(x,y) 



G T (x,y) + 



x — > x 



This completes the discussion of the quark chiral-odd channel. 



19 

-^K x - v)- 





(36) 



4 Flavour singlet parity odd sector. 

Let us now address the flavour singlet parity odd sector responsible for the evolution of axial- vector 
distribution amplitudes. For even parity there are few subtleties, which will be discussed elsewhere. 
Here we would only like to note that in the latter case a direct leading order calculation provides 
a result that suffers for the mixed channel from off-diagonal matrix elements in the unphysical 



S 



sector. Although the improved result has been found in Ref. it still remains a difficult task 
to find an appropriate representation for the dotted kernels and the two-loop G functions. 

Making use of the known non-diagonal part of the ER-BL kernel (|]), the whole NLO result in 
the axial-vector case reads 



V d)A = _y* I fy{0)A + (30 ^ _ g I y {0)A + y(0)A g D ^ G 



(37) 



where the kernels D A (x,y) and G A (x,y) are purely diagonal. Here the matrix of the LO kernels 
is given in a compact form by 



V® A (x,y) 



( 



C h 



QQ 



v(x,y) 



-2T F N f QG v a (x,y) 



\ 



^ C F GG: v a {x, y) C A [ GG v A {x, y)\ + - &5{x - y) J 



(3* 



where °- Q v = QQ v a + QQ v b and GG v A = 2 GG v a + GG v b . The general structure of the functions v 1 is 



AB„A 



v\x } y)=6(y-x) AB f\x } y)± 



x — > X 

y^y 



for 



A = B 



with 



ABja 
AB fb 



X 



v(A)-l/2 



y 



v(B)-l/2 I _J_ 



(39) 



(40) 



y—x 



The index z/(A) coincides with the index of Gegenbauer polynomials in the corresponding channel, 
i.e. u(Q) = 3/2 and u(G) = 5/2. The dotted kernels involved in the definition (|37|) can simply 
be obtained by differentiating LO results w.r.t. the index v which gives rise to the additional 
ln(x/|/)-multiplier in front of the former 



V i0)A (x,y) 



QQ 



v(x,y) 



-2T F N f QG v a (x,y) 



C F GQ v a (x,y) 



G- 



GG 



v (x,y) 



with the matrix elements 

AB 



x 



v(x,y) = 6(y - x) f(x,y)\n-± 

y 



x — > X 

y -> v 



for 



(41) 



A = B 
A^B 



(42) 



Note that for A = B the dotted kernels are defined with the "-(-"-prescription. The g function is 
given by 

~ln( 1- 

V I f 4. T* I 

(43) 



( —Ci 



g(x,y) = 9(y-x) 



y-x 



\ 







y—x 



± 



J 



X — > X 

y->y 



with (— ) + sign corresponding to (non-) diagonal elements. Note, that we have used the property 
(X - V) ln(l 
of Ref. IE7I. 



(X — X>) ln(l — |) = —(X — X>)| for the element of GQ-channel to make contact with the results 
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Next we construct the diagonal G(x,y) kernel. At first glance one would naively expect that 
one can obtain these kernels by only inserting appropriate AB f functions in the definition fl24j) , 
so that the symmetry properties of the / functions w.r.t. the weight induce then the desired 
symmetry of the AB G functions. Unfortunately, the symmetry is not sufficient for the diagonal 
form of the G(x, y) kernel. To ensure the diagonality, we have to add terms containing single logs 
and rational functions. Let us define the matrix 



G A (x,y) 



t 



\ 



2C p {C t 



2 J 



^G A (x,y)]^ 2C A T F N f Q G G A (x,y) 



C F C A G ®G A (x,y) 



C 2 



GG 



G (x,y) 



(44) 



with the following general structure 



ABri A i 



G A (x, y) = 6(y - x) ( AB H A + A AB H A ) (x, y) + 6(y - x) ( AB H A + A AB H A ) (x, y). (45) 



ABtj A 



Here analogous to the non-singlet case we set 



AB H A (x,y) = 2 



± AB f A {Li 2 {x) + \ny\nx) 
( AB f A T AB T) ( Li 2 ( 1 - 



AB f A U 2 (y) 



(46) 



x 



\n 2 y) + AB f A (Li 2 (y) - Li 2 (x) - \ny\nx) 



(47) 



where the upper (lower) sign corresponds to the A = B (A ^ B) channels. An explicit use of 
the reduction P — > V D procedure ([]) to restore the AH contributions is rather involved due to 
complexity of the integrand function. Rather we have succeeded to deduce them using different 
arguments. Since the crossed ladder diagrams have no UV divergent subgraphs the kernels AB G 
in different channels are related in a scheme independent way by supersymmetry and conformal 
covariance of M = 1 super Yang-Mills theory [Q. Employing these symmetries we restore^ all 
necessary terms in a straightforward manner to be 



A^^fx,!,) = A QQ H A (x,y) =0 



A QG H A (x,y) 
A GQ H A (x } y) 

A GG H A (x,y) 

A GG H A (x } y) 



2— In x- 2— In y, A QG H (x,y) = 2—\nx - 2^ z \ny, 

yy yy yy yy 



2 — lnx-2— \ny, A GQ H (x,y) = -2 — lnx + 2— \ny 



xx 



XX 



XX 



y 



y 



y 



x 2 1 + (x — yy 



y 



y 2y2 



„n, x x(x — x) , n x(x — x) , 
2— In - + 2^ — — — In x - 2— — — — In y, 



■r,2 



V 



y 



yy 



yy 



(48) 
(49) 

(50) 
(51) 



x x 2 1 
2 - -- + 2 — 

y y yy 



XX 



„n, x Ax + y)x, x 1 — xx, ... . , 
2— In - + 2^ ^— In 2 — In x + 6— lny. 

V x yy 2 y yy yy 



Finally, we have to extract the remaining diagonal piece D A of V A in the forward limit 



from the known DGLAP kernel P | 2q| . We take into account the underlying symmetry of the 



3 The details will be presented elsewhere. 
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singlet parton distributions to map the antiparticle contribution, i.e. z < 0, into the region z > 0. 
As expected we find from 



D A (z) = P A (z) 



LIM j-V d (V°> A + ^ - <j ® V® A + V (0)A ® </ + G A | 



(52) 



a simple convolution-type representation for ER-BL kernels which can be immediately deduced 
from the forward results for singlet QQ-channel 



^D A = C F [D F } + - C F ^ [D fi ] + - C F (c F - Cf) 



l _QQ v + 2 QQ v a 



6C F T F N f QQ v a , (53) 



where Dp, Dp are given by Eqs. (|T9| ) and fl22"|), respectively. The rest of channels is expressed as 
QG D A = ZCpTpN f { Q ^ a ® QG v a - l - QG v a } (54) 



2 C A T F N f \3 QQ v a ® Q V + [1 + 2C(2)] Q V 



OG v A 



C f Ca s - 



r 2 



GG V A- 




2 




+ 





+ - GG V 
2 



GO/ 



V + 3 



i GC V 



- - G V 

6 

5 " 2C(2) 



GQ v a 



GG v A 



Y G v a -25{x-y) 



(55) 



(56) 



'2 3 



+ G S a + 25(x - y) 
- C F T F iV / | G S a ® G V- G V + 5(x-y)}, 
where we have introduced a new kernel 



GQ v c {x, y) = 9{y - x)— {2xy - y) - { 

y [y-*y 



(57) 



These results provide us with the explicit parity odd singlet evolution kernels. 

5 Conclusions. 

In this paper, we have presented a simple method for construction of the exclusive evolution 
kernels in NLO from the knowledge of the conformal anomalies and the available two-loop splitting 
functions. The main task was, of course, the reconstruction of the diagonal part of the kernel in 
the basis of Gegenbauer polynomials. 

In the course of study we have established convolution-type formulae for the bulk of contribut- 
ing two-loop graphs with an exception of cross-ladder diagrams. The complications which arise 
in the restoration of the latter from the known forward kernels has been overcome making use of 
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M = 1 supersymmetric constraints |[22|| . The former feature suggests that by disentangling the 



topology of corresponding diagrams, it might allow for an effective and facilitated way of explicit 
calculation. One may expect that this property persists for a subset of diagrams at higher orders 
and can be used, e.g. for diagrammatical derivation of NNLO splitting functions. 

The details of the present formalism together with the flavour singlet parity even case, where 
new subtleties appear, will be discussed elsewhere. 

A.B. was supported by the Alexander von Humboldt Foundation. 
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